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SPECTRAL ESTIMATORS OF ABSORBED PHOTOSYNTHETICALLY 
ACTIVE RADIATION IN CORN CANOPIES 


ABSTRACT 


Most models of crop growth and yield require an estimate of canopy 
leaf area index or absorption of radiation; however, direct measurement 
of LAI or light absorption can be tedious and time consuning. The 
objective of this study was to develop relationships between photosyn- 
thetically active radiation (PAR) absorbed by corn ( Zea mays L. ) cano- 
pies and the spectral reflectance of the canopies . Absorption of PAR 
was measured near solar noon in corn canopies planted in a field exper- 
iment conducted at the Purdue University Agronomy Farm, West Lafayette, 
IN, on a Typic Argiaquol at densities of 50,000 and 100,000 plants/ha. 
Reflectance factor data were acquired with a Landsat MSS band radiome- 
ter. From planting to silking, the three speotrally predicted vegeta- 
tion indices examined (ratio of red to near infrared reflectance , nor- 
malized difference, and greenness) were associated with more than 95? 
of the variability in absorbed PAR. The relationship:, developed bet- 
ween absorbed PAR and the three indices were evaluated with reflectance 
factor data acquired from corn canopies planted in 1979 through 1982 
that excluded those canopies from which the equations were developed. 
Treatments included in these data were two hybrids , four planting den- 
sities (25, 50, 75, and 100 thousand plants/ha), three soil types 
(Typic Argiaquol, Udollio Ochraqualf, and Aerie Ochraqualf) and several 
planting dates. Seasonal cumulations of measured LAI and each of the 
three indices were associated with greater than 50? of the variation in 
final grain yields from the test years. Seasonal emulations of daily 
absorbed PAR were associated with up to 73? of the variation in final 
grain yields. Absorbed PAR, emulated through the growing season, is a 
better indicator of yield than emulated leaf area index. Absorbed PAR 
may be estimated reliably from spectral reflectance data of crop cano- 
pies . 
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INTRODUCTION 

Remote sensing from aerospace platforms can provide information 
about crops and soils which could be useful for modeling crop develop- 
ment and production. The feasibility of utilizing multispectral data 
from satellites to identify and measure crop area has been demonstrated 
(MacDonald and Hall, 1980), however, relatively little research has 
been conducted to develop methods of incorporating multispectral data 
into models that provide information about crop condition and yield. 
Most models of crop growth and final yield (Arkln et al., 1975; Coelho 
and Dale, 1980; Stapper and Arkin, 1979; Steven et al., 1983) require 
an estimate of a canopy's green leaf area index (LAI), absorption of 
solar radiation (SR) or, more appropriately, absorption of photosyn- 
thetically active radiation (PAR). Measurements of LAI or absorbed PAR 
are tedious and time consuming for small research plots (Daughtry and 
Hollinger, 1984) and impossible to obtain over large areas. 

Nunerous spectrally predicted vegetation indices have been pro- 
posed and used to makB quantitative estimates of LAI, phytomass, and 
percent soil cover (Asrar et al., 1984; Perry and Lautenschlager , 1984; 
Steven et al., 1983). These vegetation indices exploit differences in 
reflectance patterns of green vegetation and other materials within a 
scene. The simplest index is the ratio (Bq. 1) of near infrared (800 
to 1100 nm) and red (600 to 700 nm) reflectances. 

RATIO s IR/red (1) 

A closely related vegetation index is the normalized difference (Eq. 2) 
which is the difference in the ln:’r?u- d and red reflectances divided by 
the sum of the two reflectances . 

ND s (IR - red) / (IR ♦ red) (2) 

The third vegetation index is the *p?eer.ness index (GI) which is a 
linear transformation that accentuates the presence of green vegeta- 
tion (Rice et al., 1980). The greenness index (GI) is computed as: 

GI = -0.4894 RF 1 - 0.6125 RF 2 + 0.1729 RF^ + 0.5954 RF^ (3) 

where RF. through RFj, are the reflectanco factors in bands 500 to 600, 
600 to 7o0 , 700 to 800 and 800 to 1100 nm, respectively. While GI is 
sensitive to green vegetation it is relatively insensitive to changes 
in the amount of shaded area within row crop canopies that result from 
the diurnal changes in solar azimuth and zenith angles (Kollenlcark et 
al., 1982). 

Leaf area index and the proportion of solar radiation intercepted 
by corn canople' may be estimated with GI (Daughtry et al., 1983). 
Similarly, LAI aiiu the proportion of PAR absorbed by wheat can be esti- 
mated with ND or GI (Asrar et al., 1984; Hatfield et al., 1984). 
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The seasonal duration of leaf area Is frequently a more Important 
Indicator of grain yields than maxlmm LAI produced. Correlations of 
grain yields with speotral or agronomic data aoqulred on a single date 
during the growing season of a crop may be spurious and must be used 
with caution (Daughtry et al., 1983). Vegetation Indices accunulated 
for various portions of the growing season have been found to be asso- 
ciated with significant proportions of the variance in grain yields of 
wheat (Pinter et al., 1981; Tuciosr et al 1980), corn (Daughtry et 
al., 1983; Walburg et al., 1982) and sugar beets (Steven et al., 1983)* 

The daily Interpolation of spectral transformations , or spectrally 
predicted agronomic variables, through a crop's growing season have 
been related to grain yield. The normalized difference computed with 
the 650 to 700 run (red) and 780 to 820 nm (near IR) wavebands and 
interpolated for portions of the growing season was associated with 64$ 
of the variance in grain yield of wheat (Tucker et al., 1980). 

Daily dry matter production may be described as a function of the 
total incident solar radiation, the proportion of solar radiation 
intercepted, and the efficiency of conversion of solar energy to phyto- 
mass (Steven et al., 1983). When the proportion of intercepted SR was 
estimated with multispectral data and the predicted dally dry matter 
production was cumulated over the growing season, Steven et al. (1983) 
were able to predict dry matter of sugar beets at harvest within 6$. 

Our study was conducted to develop relationships between absorbed 
PAR and several vegetation indices and to examine the relationship bet- 
ween grain yield and spectrally estimated absorbed PAR. 
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MATERIALS AND METHODS 

This -tudy consisted of experiments conducted over four years at 
the Purdue University Agronomy Farm, West Lafayette, IN (4oP 28' N, 87° 
00' W). Two replicates of an adapted corn ( Zea mays L. ) hybrid were 
planted on several dates at several plant densities (Table 1). The two 
soil types in Experiments I and II had distinctly different spectral 
reflectance characteristics in the red and near IR wavebands. The corn 
was planted in north-south rows with 76 cm spacing between rows. Soil 
analyses were conducted and N, P, and K applied to maintain high levels 
of fertility. Preemergence herbicides were applied for weed oontrol. 

Daily meteorological data were recorded at the cooperative 

National Weather Service station (Wesc Lafayette 6NW) that was within 
200 m to the east of the plots . Incoming solar radiation was measured 
with an Eppley Precision Spectral Pyranometer and was recorded as total 
MJ m -2 day" . The photosynthetically active portion of the daily SR was 
assuned to be 0.5 (McCree, 1966; Szeicz, 1974). Daily incident photo- 
synthetic photon flux density (PPFD) was oomputed as: 

PPFD = 0.5 (SR) 4. 6X 1 0“ 6 (4) 

where SR is the daily incident solar radiation end 4.6X10"^ is the con- 
version factor (moles/J) determined by McCree (1972). 


Canopy Characterization 

Agronomic variables were usually measured at weekly intervals and 
included leaf area index and development stage (Ritchie and Hanway, 
1982). Green leaf area was measured with a LI-C0R model LI-3000 
(LI-3100 in 1982) area meter from a subsample of the plants harvested 
each sample date. Three plants were harvested on each sample date in 
1979, four in 1980 and five in 1981 and 1982. Green leaf area index 
was oomputed for each plot as the ratio o green leaf area to soil 
area. After physiological maturity, as indicated by black layer forma- 
tion, grain was harvested. Harvested areas included the middle rows of 
each plot. Grain yield was computed on an area basis (Mg/ha) and 
adjusted to 15.5% moisture. 
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Table 1. Summary of the experimental conditions for the four years that com- 
prised the study of spectrally predicted estimators of APAR'. 


Experi- 

ment 

Number 

Year 

Planting 

Dates 

Planting 

Densities 

Hybrid 

Soil Types 




3 

10 plants /ha 



I 

1979 

2, 16, 30 May 

25, 50, 75 

Beck 65X 

Typic Argifquol 
Udollic Ochraqualf 

II 

1980 

7, 22 Mr.y 
1 1 June t 

25, 50, 75 

Beck 65X 

Typic Argiaquol 
Aerie Ochraqualf 

III 

1981 

8, 29 May 
11, 29 June 

25, 50, 75 

Beck 65X 

Typic Argiaquol 

IV 

1982 

14 May 
8, 24 June 

25, 50, 75, 
100 

Adlers 3CK 

Typic Argiaquol 

V 

1982 

14 May 
24 June 

50, 100 

Adlers 3CK 

Typic Argiaquol 


+ Additional plantings of the 50,000 plants/ha density were planted on 16 and 
29 May, 18 June and 3 July 1980. 
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Measurements of Absorbed PAR 


Photosynthetic photon f ux densities (PPFD) were measured under 
clear sky conditions (cloud cover less than 10? with no olouds within 
1C P of sun) with a line quantun sensor (LI-COR 191SD;. The sensor has 
a cosine corrected response to spatially average the Incident PPFD over 
its 100. 0 X 1.27 cm rectangular surface. The sensor was modified with 
the addition of a handle and two bubble levels (one on top and the 
other on the bottom of sensor). A switch on the handle of the sensor 
allowed the observer to trigger automatic data acquisition by a data 
logger (Omnidata Folycorder Model 516). The time of each measurement 
(hour, min, s) was also automatically recorded. The sensor was always 
leveled and positioned such that no shadows from the handle or observer 
influenced measurements. Care was also taken to minimize possible 
reflectance from the observer. Transmitted PAR (TPAR), reflected PAR 
from the canopy and soil surface (RPAR^), and reflected PAR from the 
soil surface under the canopy (RPAR ) were computed as proportions ol 
the incident PAR (PAR q )s 


TPAR' = (TPAR/PAR q ) , 
RPAR' os s (RPAR Cfl /PAR 0 ) , 
RPAR' s = (RPAR a /PAR 0 ) , 


(5) 

( 6 ) 
(7) 


and will be cited as such throughout the following discussion unless 
otherwise noted. 


PAR and TPAR' were measured in Experiment V at weekly (when pos- 
sible) intervals throughout the growing season under clear sky condi- 
tions within 0.5 hr of solar noon. TPAR' was measured at three sites 
per plot. At each site the sensor was positioned perpen icular to the 
row direction ar.d centered on the row. Each set of TPAR measurements 
consisted of four individual measurements as the sensor was incremented 
at nearly equal intervals between plants within a row. TPAR' was com- 
puted as the mean of four individual measurements and PAR (Eq. 5). 
PAR q was measured either above or outside of the canopy within 20 s of 
themeasured TPAR. TPAR' measured with the 100 cm sensor length varied 
from that of the ideal sensor length (Warren Wilson, 1981) of 76 cm by 
less than 71 when less than 8 leaves were emerged (stage V8) . After 
this stage leaves overlapped between rows and few significant differ- 
ences were detected. 


RPAR fl and RPAR^ varied less than 6? during the growing season. 
RPAR was measured with the sensor Inverted and leveled 35 cm above 
the mean height of the canopy at two sites per plot. One site centered 
the sensor over and perpendicular to the crop row. The second site cen- 
tered the sensor over and perpendicular to mid rcw. RPAR' a is the por- 
tion of the PAR transmitted through the canopy that is reflected from 
the soil surface. Direct measurement of RPAR was net feasible as 
placement of the line sensor under the canopy at a neight of 2.0 cm 


'( 


above the soil surface (Hlpps et al., 1 983) resulted in a shadow oast 

on the soil surfaoe by the sensor. RPAR' was estimated as: 

9 

RPAR* 8 = (RF S )(TPAR') (8) 

where, RF v is the refleotanoe faotor of bare soil. RF was measured 
with the sonsor inverted and leveled 35 cm above the sol? surfaoe at a 
site adjacent to the plots. RF of the dry soil (Typic Argiaquoll) was 
measured as 9.3 + 0.551 (n*270) . RPAR' (Eq. 8) ranged from 0.1 of PAR 
when no oanopy cover was present to less than 0.01 of PAR under a ful? 
canopy oover. Absorbed PAR v ”q . 9) was oomputedt 

A PAR' s 1.0 + RPAR* - TPAR' - RPAR' . . (9) 

9 09 

Speotral Measurements 

Radianoe measurements, used to determine reflectance factors (RF), 
were acquired with a four band radiometer (Exotech model 100), that 
simulates the Landsat MSS bands (MSS4, 500 to 600 nm; MSS5, 600 to 700 
nm; MSS6, 700 to 800 nm; MSS7, 800 to 1100 nm). Measur vnents were made 
throughout eaoh growing season at approximately weekly intervals. 
Biehl and Robinson (1983) describe the conditions and procedures for 
obtaining the RF data. The radiometer has a Iff field of view. 

The radiometers were attached to a boom mounted on a pickup truck 
and elevated 7.6m (5.2 m in Experiment I) above the soil surface. 
Data were taken only when there were no clouds in the vicinity of the 
sun and when the solar elevation was at least 4^. Measurements were 
made after the instruments were leveled for a nadir view angle. Each 
set of measurements consisted of two measurements, one centered over a 
row and one centered between rows of the canopy to better estimate the 
overall canopy response (Daughtry et al., 1982). Three seta of measure- 
ments were made for each of the plots included in Experiment V. 

Analysis of Data 

APAR' measured in Experiment V was regressed as a function of 
green LAI and each of the three speotral transformations . The rela- 
tionships developed were tested on the data collected in Experiments I 
through IV (Table 1) The daily amount of APAR (moles m~ 2 ) was estimated 
as: 


APAR = APAR' PPFD (10) 

where APAR' is the predicted portion of the daily incident PAR (PPFD) 
that was absorbed by the canopy. 

The green LAI and vegetation indices , their respective predicted 
values of APAR' and estimated values of APAR (moles in ) were linearly 
interpolated and emulated from planting to maturity for each field 
plot Included in the test data. One plot in Experiment I and eight in 
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Experiment II lodged severely end were exoluded. Representative values 
of refleotanoe from bare soil and senesoed vegetation were used respec- 
tively to begin and end seasonal Interpolations of APAR* . Final grain 
yield and the results of the seasonal emulation of agronomio and vege- 
tation indices and their estimates of APAR' were averaged over the 
replioates. The variance in final grain yield associated with the 
emulated agronomio and vegetation indices was examined. 
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RESULTS AND DISCUSSION 


Relation of APAR' to LAI and Spectral Variables 

Two distinct relationships between APAR' and LAI were detected 
(Fig. 1). APAR' increased as a function of green LAI from planting to 
a maximum at anthesis or silking (stage Ri) then decreased at a diffe- 
rent rate to maturity. Beer's Law described the relationship between 
APAR' and LAI from planting to silking (Asrar et al., 1984; Hipps et 
al. f 1983; and Norman, 1980). Similar relationships between APAR' and 
LAI were reported by Hatfield et al. (1984) for wheat canopies. Maximum 
APAR' for wheat canopies also occurred near anthesis. The different 
relationships before ana after silking are due to the absorption of PAR 
by nongreen vegetation after silking as senescence occurs. Relation- 
ships developed after silking between APAR' and LAI (or the vegetation 
indices) overestimate the APAR' of green vegetation and hence the uti- 
lization of incident PAR by the conchy. Even at maturity, when no green 
leaves were present, greater than 60} of the incident PAH was absorbed 
by nongreen leaves, stems, and other plant parts (Fig. 1). 

Quadratic equations (Table 2) adequately described the relation- 
ship between APAR' and the spectral variables RATIO (Fig. 21, ND (Fig. 
3), and GI (Fig. 4). Only data obtained between the planting and silk- 
ing stages of crop development were included in the development of 
these equations. The equations developed in Experiment V were tested 
with the data of Experiments I through IV. APAR' predicted from mea- 
sured LAI was highly correlated with APAR' predicted with the three 
vegetation indices (Table 3). 

A canopy with north-south row orientation has a minimum APAR' at 
solar noon (Hipps et al., 1983). Thus daily estimates of the quantity 
of absorbed PAR (moles/m2) based on the proportion of absorbed PAR 
(APAR') measured at solar noon probably under-predict the energy avai- 
lable to the crop before a full canopy has developed. When LAI of 
wheat is high (i.e., > 4.0), the absorption of PAR remains high 
throughout the day and is largely independent of solar angle. In wheat 
canopies with low LAI, proportionately more PAR was absorbed on cloudy 
days than on clear days (Hipps et al.,1983). Nevertheless Hipps et 
al.(1983) reported that a single equation adequately described absorbed 
PAR as a function of LAI. 



0 1 2 3 4 5 6 7 8 


LAI 

Figure 1. Relation between absorbed PAR and leaf area index for growth 
(planting to silking) and senescent (silking to maturity ) periods of 
corn (Experiment V). The solid line indicates the predicted values of 
APAR' for the interval from planting to silking. 
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RATIO 


Figure 2. Relation between absoroed PAR and RATIO for growth (planting 
to silking) and senescent (silking to maturity) periods of corn devel- 
opment (Experiment V). The solid line indicates the predicted values of 
APAR' for the interval from planting to silking. 

!> 


ZJ 



Figure 3« Relation between absorbed PAR and normalized difference (ND) 
for growth (planting to silking) and senescent (silking to maturity) 
periods of corn development (Experiment V). The solid line indicates 
the predicted values of APAR' for the interval from planting tc silk- 
ing. 



Figure 4. Rel-.' .on between absorbed PAR and greenness index (GI) for 
growth (planting to silking) and senescent (silking to maturity) per- 
iods of corn development (Experiment V). The solid line indicates the 
predicted values of APAR' for the Interval from planting to silking. 


12 


Thble 2. Quadratic regression equations that were developed from data of Exper- 
iment V to predict APAR', from planting to silking, with spectrally predicted 
vegetation indices (n=40). 


estimator 


coefficients 

RMSE 

F 

R 2 


b 

b 

b 

(*) 



GI 

-C. 20 + 

0.058 

-0.0007 

5.3 

485.2 

0.96 

ND 

0.6 

-2.2 

2.9 

4.7 

613*7 

0.97 

RATIO 

-0.06 

0.102 

-0.0026 

4.4 

709.5 

0.97 


+ All coefficient were significant at the 0.01 level of probability, except b 
of the RATIO quation, which was significant at 0.1. 


Table 3. Correlation coefficients that resulted from the prediction of APAR' 

for the test years (Experiments I through IV) with the equations developed 
from the data of Experiment V (n=56l). 


predictors of APAR' 



LAI 

GI 

ND 

RATIO 

LAI 

1.00 + 




GI 

0.96 

1.00 



ND 

0.92 

0.91 

1.00 


RATIO 

0.94 

0.94 

0.99 

1.00 


t 


All correlations were significant at the 0.01 level of probability 
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Relation of Yield to Cumulative LAI and APAR 

The seasonal duration of leaf area Is a better Indicator of grain 
yields than maximum LAI produced (Daughtry et al . , 1983)* Seasonal 
emulations of dally LAI or any of the three vegetation Indices were 
associated with greater than 50} of the variation In oorn grain yields 
(Table 4). LAI of a crop represents only the amount of photosynthetic 
tissue present In the canopy and does not account for productivity. 
The use of only daily incident PAR (Fig. 5a) or APAR' (Fig. 5b) in a 
model of crop productivity clearly misrepresents the actual amount of 
PAR absorbed (Fig. 5c). emulated daily APAR was associated with over 
66} of the variation in final grain yields measured from four years of 
test data (Table 4). The emulated daily APAR ( ZAPAR) , computed with 
APAR' predicted from the normalized difference (Eq. 2) and daily inci- 
dent PAR, was associated with 73} of the variation in observed grain 
yields (Fig. 6). A portion of the unexplained variation in grain 
yields is due to the effects of meteorological variables not included 
in this "odel on plant growth and development . The various planting 
date and plant density treatments would also be associated with a por- 
tion of the mexplained variation in the grain yields . 

In summary, the intent of this study was to develop and examine 
the performance of several spectral estimators of one specific variable 
often included in models of crop yield (APAR 1 ). The relationships bet- 
ween final grain yield and APAR indicate that all three of the vegeta- 
tion indices that estimate APAR' from canopy spectral reflectance pro- 
vided similar information about the portion of incident PAR that was 
absorbed by crop canopies. Spectral estimators of APAR, interpolated 
on a daily basis and emulated from planting to maturity, were associ- 
ated with as much of the variation in final grain yield as APAR esti- 
mated with measured LAI. These results suggest that APAR may be esti- 
mated from canopy spectral reflectance for large areas where direct 
measurements of LAI would be prohibitive. Thus models of crop yields 
which require estimates of absorbed PAR may be implemented and evalu- 
ated . 


APAR (moles/m 1 ) APAR 1 INCIDENT RARlmoles/m*) 



CALENDAR DAY 


270 


. 5a. Daily incident PAR for 
the 1982 growing season. 




Spectrally predicted (green- 
ness index) portion of inci- 
dent PAR absorbed by corn 
(APAR') planted on I 1 ! May 
1982 at 100,000 plants/ha. 
The development stage of the 
corn is indicated for 
selected observation dates. 


\ 


Daily amount of PAR absorbed 
by a corn canopy (product of 
daily incident PAR and APAR') 


290 
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Table 4. Results of linear regression of corn grain yields of Experiments I 
through IV on seasonally cumulated values of agronomic and vegetation indices 
and their respective estimates of canopy absorption of PAR (n=79) . 


variable 

RMSE 

F 

r2 


Mg/ha 



LAI 

2.3 

98.4 

0. 56 

GI 

2.3 

100.4 

0.56 

ND 

2.2 

105.6 

0.58 

RATIO 

2.3 

101.1 

0.57 

zap “lai 

1.8 

193-3 

0.72 

EAPM C i 

2.0 

151.3 

0.66 

:apa» nd 

1.8 

209.2 

0.73 

EAPAR RATIO 

2.0 

146.4 

0. 66 


20 

18 

16 

14 

*o 

^ 12 

o» 

I .0 

Q 

_l 

UJ 8 


r*=0.73 
RMSE = I.8 
n= 79 
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Figure 6. Grain yields of corn from Experiments I through IV as a function 
of absorbed photosynthetic photon flux density (ZAPAR). 
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